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It has been predicted that in the presence of a sufficiently high-dissipative environment transport 
in a small tunnel junction can become extremely regular, giving rise to the phenomenon of single- 
electron tunneling oscillations. Recent progress in detection of high-frequency current fluctuations 
and the interest in single-electron sources motivate further investigations on the expected accuracy 
of the charge oscillations as a function of the impedance of the environment. In this paper we study 
theoretically the charge-fluctuation spectrum at finite frequency for the system at hand, and inves- 
tigate its behavior as a function of the external impedance. The evolution and the disappearance 
of the single-electron oscillations peak is described by analytical and numerical methods. 



I. INTRODUCTION 

Coulomb blockade of voltage-biased tunnel junctions 
in the presence of a dissipative environment is a well- 
understood phenomenon of quantum transport (for a re- 
view see for instance Ref. [l|). It has been shown that if 
the impedance of the environment R is smaller than the 
tunnel junction resistance Rt, but larger than the quan- 
tum of resistance Rq = 2nh/c w 25.8/cfi (H being the 
reduced Planck constant and e the electron charge), the 
current is suppressed for bias voltage V smaller than the 
Coulomb gap e/2C, associated with the capacitance C of 
the junction (C is typically in the fF range). This effect 
can be visible for temperatures lower than the Coulomb 
charging energy e 2 /2C and for R t 3> Rq, which assures 
suppression of cotunncling. In practice realizing a high- 
impedance environment with a flat frequency response 
till frequencies of the order of e 2 /2CH is a challenging 
experimental problem, as discussed in detail in Ref. 0. 
Most experimental observations of Coulomb blockade 
phenomena in single tunnel junctions are actually done 
in an intermediate-impedance situation (R ~ Rq), which 
leads at least to clearly non-linear characteristics of the 
junction; see for example Refs. @, 0|. Nevertheless, it is 
feasible to realize impedances of the order of some hun- 
dreds of kO; for instance in the recent Ref. @ a resis- 
tance R ~ 0.4M17 has been realized. For these values not 
only the suppression of the current at low-bias voltage 
should be visible, but also an appealing effect predicted 
m the 1980s mm-- the single-electron tunneling oscilla- 
tions (SETOs). 

The idea behind this effect can be understood in the 
simplest way in the limit R, V — > oo with V/R = lb, so 
that the tunnel junction is current biased. The current 
slowly charges the capacitance. When V = Q/C reaches 
the threshold e/2C one electron can cross the junction. 
If <C e/R t C this will happen just after V has reached 
the threshold. The charge on the capacitance after the 
tunneling event will be Q ~ — e/2, and Q will start to 
slowly increase again. A time ~ e/Ib is needed before 
a new electron can cross the junction and the sequence 
can start again. The voltage at the junction will thus be 



periodically modulated at a tunable frequency ~ h/e. 

Observation of this phenomenon is difficult. The im- 
plementation of the required strong-impedance environ- 
ment has been realized by different authors using on- 
chip resistors [Tfj| - [l6j . An alternative approach has also 
been tried by designing the environment with tunnel- 
junction arrays [l?], EH, exploiting the fact that a large 
number of arrays reduces the stochastic nature of the 
current [l9j]. Recently, observation of soliton-like single- 
electron oscillations with this method has been reported 
[20l | . A similar phenomenon for superconducting Joseph- 
son junctions has been predicted (Bloch oscillations) [21[ 
and investigated by many authors [22rl25| . Arrays of dc 
SQUIDs have also been exploited in this case to build 
up the proper environment to obtain Coulomb block- 
ade of Cooper pairs [2(j| - |28j . Reports on the observation 
of Shapiro-step-likc structures in microwave-irradiated 
junctions constitute the present state of the art for the 
experimental probe of this effect [HI, [H, H^] . 

Progress in the detection of high-frequency current 
fluctuations f30l - l33l ] can open new possibilities of obser- 
vation of this phenomenon and of the crossover region, 
where the oscillations are not completely established. At 
the same time the possibility of generating a periodic and 
frequency-tunable electric signal without any oscillating 
source is an interesting opportunity and could have ap- 
plications, for instance, as a motion actuator in nanomc- 
chanical systems or a controlled single-electron source. 

It should be mentioned that Coulomb blockade is in- 
stead easily observed in double tunnel junctions even in 
absence of an environment. This happens since the sec- 
ond tunnel junction plays the role of the large-impedance 
environment suppressing quantum charge fluctuations 
and preventing tunneling if the voltage is below a thresh- 
old. Note however that this configuration (without envi- 
ronment) will not give rise to single-electron oscillations, 
but to the well known sequential transport regime. Elec- 
trons hop stochastically through the first and then the 
second junctions: The advantage of an ohmic environ- 
ment is that it can generate (at least in principle) a stable 
current source. 

In this paper we study how accurate the SETOs can be 
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FIG. 1. Circuit scheme of the system considered: a tunnel 
junction with capacitance C and tunneling resistance Rt, bi- 
ased by the constant current and shunted by a resistance 
R s in parallel (right), which is equivalent to a junction bi- 
ased by a voltage Vb and with a resistance R in series (left) 
provided that Vb = hR and R = R s - 



as a function of the impedance and of the bias conditions. 
In order to do this we will study the charge-fluctuation 
spectrum at the junction capacitance (or cquivalently 
the current-fluctuation spectrum through the resistance 
load), and in particular the width of the peak at the 
frequency ~ i&/c. Notwithstanding the relatively large 
number of papers on this subject, some specifically ad- 
dressing the dynamics of arbitrarily biased mesoscopic 
tunnel junctions with an analytical approach j33 - l3o| . a 
consistent calculation of these quantities is not available. 
It can be useful in order to evaluate the expected effect 
in view of measuring current fluctuations in this kind of 
device. In this paper we show that the width of the peak 
scales as the inverse of the impedance for large R/Rt- 
The peak remains observable till values of R/Rt of the 
order of 5. 

The plan of the paper is the following. In Sec. [TT| we 
present the model describing the transport through the 
junction. In Sec.|TTT]we discuss the different regimes that 
the system undergoes by varying the bias voltage and the 
environment resistance. In particular we obtain an an- 
alytical expression for the the I-V characteristic in the 
SETOs regime. In Sec. IIVI we calculate analytically the 
charge spectrum. The results are discussed and com- 
pared with numerical Monte Carlo simulations in Sec.fVl 
Section PVTl gives our conclusions. 



II. THE SYSTEM AND THE MODEL 

Let us consider a tunnel junction with tunneling re- 
sistance R t ^> Rq and associated capacitance C. The 
circuit is voltage biased (at voltage Vb) at zero tempera- 
ture (T = 0) in the presence of a resistor of resistance R 
in series with the junction (see Fig. (UJ, left side). This 
circuit is equivalent to one with a current source I b and 
a shunt resistor R s in parallel to the junction, provided 
I b = V b /R and R = R s (see Fig. (JT), right side). We 
will thus use the parallel configuration to describe the 
device in analogy with the previous literature [7j- It is 
clear that the results can be readily converted to the volt- 
age bias case. In particular note that the limit R s — > oo 
(specifically R s ^> R t ) describes the ideal current source. 



Since we are interested in studying SETOs we need 
R s > Rt S> Rq- We thus assume from the outset that 
R s 3> Rq, which allows to neglect quantum fluctuations 
and treat the charge degrees of freedom classically 0]. 
We also assume that the environment has a flat frequency 
response Z(oS) = R s up to frequencies hu pa e 2 /2C. This 
hypothesis, though not easy to fulfill in practice, is the 
common assumption in the literature about this problem, 
and allows a simpler and more transparent approach. In 
this regime transport through the junction is described 
by the theory currently known as orthodox Coulomb 
blockade theory. Specifically the electron-tunneling rate 
depends on the voltage at the junction (Vj) as follows: 

T(Vj) = 9(Vj - e/2C)(Vj - e/2C)/(eR t ) , (1) 

where 9 is the Heavisidc function. [At finite temperature 
the function 9(V) is substituted by l/( e - cV ^ BT - 1).] 
If (Rq <c) R s <C Rt the standard picture of Coulomb 
blockade applies to the degrees of freedom of the en- 
vironment: They have the time to relax to thermal 
equilibrium between two electron-tunneling events. The 
current-voltage characteristic in this case is then given 
by Ij{Vj) = (Vj - e/2C)9(Vj - e/2C), exposing a clear 
Coulomb gap for the current Ij through the junction. 
But if R s ~ Rt or larger (for moderate values of the bias 
voltage) the resistive environment (described for instance 
by a large collection of bosonic modes) reaches thermal 
equilibrium, but not the charge on the capacitance Q(t), 
which needs a time r s = R S C to relax to its stationary 
state. Formula ([TJ still holds, but with a time-dependent 
voltage V(t) = Q(t)/C. The time dependence of the 
charge is given by the solution of the differential equa- 
tion: 



Q = -Q/RsC + I b 



(2) 



which for an initial condition Qo at t = reads 

Qf(Qo,t) = {Qo - IbT a )e-*' T ' + I b r a . (3) 

The stochastic problem is then completely formulated 
and in the remainder of the paper we discuss the behavior 
of the current and of the charge as a function of the 
two relevant dimensionlcss parameters of the problem: 
p = R s /R t and « = (J 6 — J t h)/^th, with I t h = e/2r s 
the threshold for the current to start flowing through the 
tunnel junction. 

The current through the junction has already been cal- 
culated numerically in the very early literature @. We 
obtained the same results by Monte Carlo simulation 
and, for convenience, we reproduce the curves in Fig. 
The limit of infinite R s , or ideal current source, was dis- 
cussed in details in Rcfs. 0,113] • There it was shown in 
particular that in this limit the system is in the SETOs 
regime with frequency I b /e with an averaged voltage at 
the junction given by (see 0]): 
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FIG. 2. (Color online) Average current through the junction 
versus average voltage for different values of the ratio p = 
Rs/Rt- The curves evolve from the standard Coulomb block- 
ade suppression (rightmost line corresponding to p — 0.1) to 
a square root behavior (leftmost line p = 500). The dotted 
line is the large- p limit of Eq. Q. 



In the next section we will discuss the behavior of the sys- 
tem for the intermediate regimes appearing when R s /R t 
is not infinite, deriving in particular new analytical ex- 
pressions for the SETOs frequency and I-V characteris- 
tics. 



III. REGIMES OF CURRENT TRANSPORT 

In this section we study the evolution of the current 
through the junction as a function of the current bias for 
different values of the external resistance. The most in- 
teresting case is when R s /R t is very large, we thus plot 
in Fig. [3] the current on a logarithmic scale for the ex- 
treme value of Rs/Rt = 5 x 10 2 . Fig. ((3)) will be used 
as a "map" for the rest of the section, where we will dis- 
cuss how the junction evolves through the four different 
regimes of transport indicated in the figure with roman 
numerals from (I) to (IV) . Fig. ((4| shows the behavior of 
Q(t) in the different regimes. 



A. Non SETOs regimes 

Let us begin with the region indicated with (I) in 
Fig. For I b < 7 th the whole current flows through 
the shunting resistance and the voltage at the junction 



(V 3 )i = Rsh 



(5) 



remains below the threshold of the Coulomb blockade: 
e/2C. Note that this first branch of the I-V characteris- 
tic in Fig. © is flattened on the Ij = value and it is 
thus not visible. 

For It > 7th transport through the junction becomes 
possible. For 1^ — J t h — > + one can identify a Poissonian 
regime of transport (region (II) in Fig. ((31)), where the 
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FIG. 3. (Color online) Average voltage through the junction 
versus bias current in logarithmic scale (p = 5 x 10 2 ): Four 
different transport regimes can be outlined. The solid line 
gives the results of the numerical simulations, while dashed 
lines represent the analytical curves given in Sec. IIIII 
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FIG. 4. (Color online) The time behavior of the charge in 
the different regimes. Extreme values of the parameters p and 
K have been chosen to best outline the differences. 



time between tunneling events fluctuates strongly. This 
is due to the fact that r(V) as given by Eq. ((TJ) vanishes 
linearly near the threshold, and for very small 7;, — 7th 
the charge has always the time (typically t s ) to reach 
the saturation value Q s = ItT s . The typical inverse time 
between two tunnel events is: 



1/7-eff = r(Q B ) = {R s /R t ){I b - 7 th )/e « l/r s 



(G) 



The last inequality sets also the region of existence of 
the regime (II) i.e. , < k <C 1/p. The average Vj can 
be readily evaluated by averaging the oscillations of the 
charge Qit) = Q s — ee"'' Ts on the average time between 
two tunneling events r e a + T s ~ T e s. This gives: 

(Vj) n = R s [I b (l - R s /Rt) + hhRs/Rt] (7) 
and the curve is shown dashed in Fig. ([3]). Note that 
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the slope changes sign at R s /Rt = 1- We will see that 
for R s <C Rt this region joins continuously region (IV) 
without the appearance of region (III). 

It is thus convenient to discuss now the region (IV) 
defined as the limit of large if,. In this limit the junction 
has a I-V characteristic of a normal resistor shifted by 
the Coulomb gap. The average Vj reads then: 



{V 3 )iv = 



R s Rt 
Rs + Rt 



2R t C 



(8) 



This expression holds for (Vj) 3> e/2C, i.e. for l&/ith 3> 
(Rt + R s )/R s - For R s ^> R t it is then clear that a large 
region defined by the condition 



(9) 



exists between region (II) and region (IV). This is the 
SETOs region, (III) in Fig. ([3]), which will be discussed 
below. On the other side, for R s <C Rt, one sees that 
region (II) and region (IV) overlap at re w 1. Actually it 
is straightforward to check that Eq. ([8]) expanded to first 
order in R s /Rt coincides with Eq. ([7|). 



B. SETOs regime 

Let us now discuss the single-electron oscillations 
regime, defined as region (III) in Fig. ([3]). This region 
is present only if R s R t and is characterized by nearly 
periodic electron tunneling events, since the time between 
two events is dominated by the deterministic charging 
time of the capacitance. This time is typically of the 
order of , defined as the time needed to charge the ca- 
pacitance from Q = — e/2 to Q = c/2: 
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(10) 



The electrons hop just after the threshold voltage has 
been reached. 

A general statistical theoretical framework is presented 
in Ref. Q (and recalled in the Appendix), but its ana- 
lytical solution is given there only in the ideal current 
source limit (R s /Rt — > oo). Actually in the case of the 
SETOs the approach is simplified and further progress 
is possible. In order to obtain these results and for the 
calculation of the correlation function of the next section 
it is convenient to introduce a few concepts. In Fig. §5§ 
the typical time dependence of the charge Q(t) in the 
SETOs regime is shown. We can associate a number n 
to each hopping event and define t n and t„ = t n + Sr n as 
the instant of time when Q(t) = e/2 and when the hop- 
ping event takes places, respectively (see Fig. (J5J). These 
quantities fluctuate randomly, but a correlation between 
t n and t n -i exists. Inversion of Eq. ([3]) gives the time 
needed to reach the border of the Coulomb blockade re- 
gion starting from a charge Qq: 



S(Qo) = -t s In 
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FIG. 5. (Color online) Details of the tunneling process: t n 
is the time at which the charge reaches the blockade region 
border after the (n — l)th tunneling event and St„ is the time 
it stays outside the border before the nth event. 



The following relation between successive times t n holds: 



1 1, 



(12) 



with 



3(8r) = 5T + E[Q f l-,6T 



= St — T s In 
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(13) 



Let us now introduce the probability P n (t) that n elec- 
trons have tunneled through the junction at the time t. 
Within the SETOs region, this quantity is different from 
zero only in a small time region of the order of £+, and, in 
particular, the above mentioned condition on the typical 
hopping time implies that P n (t) will reach 1 and then 
vanish in a time much shorter than i*, just after Q(t) 
crosses the threshold e/2. The rate equation for t > t n 
[Q(t n ) — e/2] takes the simple form: 

rlP 

-^(t) = -T(Q f (e/2,t)/C)P n (t) (14) 
with the initial condition P n (t n ) = 1. The solution reads: 

P n {t) = e ciit ° \ ra ) . (15) 



For short times (t — t n <C t s ) it has a Gaussian form 

P n (t) w exp {-(t - t n ) 2 K P /{AT 2 a )} (16) 

with a decay time scale ~ r s / \/kp <^ in region (III). 
The Gaussian form will thus be used in the following for 
the analytical calculations. 

From the knowledge of P n (t) it is possible to obtain 
the probability density that a hopping event takes places 
at time t: V(t) = —dP n /dt, for t > t n . This allows us 
to calculate the average delay time (St) for an electron 
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to hop after the threshold e/2 has been crossed by the 
charge Q(t): 



0.06 



(St) 



dt t V(t) . 



(17) 



In particular when P n can be approximated by the Gaus- 
sian (fTO)) one obtains (St) /t s = ^JtF/Jkp) with (St) <C i*. 
To obtain the period of the SETOs one has to average 
the nonlinear expression (|13|) : T = ($}, which again in 
the SETOs region simplifies to 



T = i* + 2 (St)/(2 + k), 



(18) 



Let us now come back to the probability. Conservation 
of the probability gives that P n+ i(t) = 1 — P n (t). Since 
in this approximation the behavior is quasiperiodic, the 
charge on the capacitor for the n+1 electron is on average 
Qf(e/2, t — t n — T). (A more precise discussion on the 
validity of this last average can be found in the Appendix, 
where the problem is analyzed more rigorously.) The 
average charge can then be computed by averaging over 
a period as follows: 



(Q) = 



t n +T 



Q/(|,i-t„)P„(t)+ (19) 
Q/(|,i-i„-T)P„+i(t) . 



In the limit of p 3> 1 only the gaussian part of the prob- 
ability is relevant, and the integral gives: 



(Q) =Q S - ( Qs - g 
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(20) 
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Eq. dSHJ) leads to (Vj)ui = C(Q). With little loss in the 
accuracy Eq. (|20|) can be simplified to the form: 
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(21) 



The analytical expressions (|20l) . (|21j) . and (|A8[) obtained 
in the appendix, are compared to the Monte Carlo re- 
sults in Fig. ©. These expressions describe the current 
with good accuracy, and in particular they all capture 
the presence of a minimum in the voltage Vj. This min- 
imum signals the crossover region between two different 
kinds of SETOs. We indicate them in Fig. © as III.l 
and III. 2. The latter appears for 1 <C K <C p. In this 
case, and in the extreme limit p — > oo the SETOs period 
becomes T = e/h, he. corresponds exactly to the time 
needed to the ideal current source to furnish a charge e. 
The saturation value for the charge (Q s ) in this regime is 
much larger than e/2, implying that only the linear part 
of the exponential in Eq. is explored. This is impor- 
tant since the small fluctuations in the hopping times do 
not affect the evolution equation for the next electron. 
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FIG. 6. (Color online) Comparison between the differ- 
ent approximations for the SETOs "nose" of the current- 
voltage characteristics: Monte Carlo data (circles) are shown 
against analytical calculations from Eq. (|A8f) (solid), Eq. (|20|) 
(dashed), and Eq. (f2"T) (dot-dashed). 
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FIG. 7. (Color online) Scheme of the boundaries of the 
transport regimes as a function of the relevant parameters: 
The tunneling/shunt resistances ratio 1/p = Rt/R s and the 
relative distance of the current bias from the threshold k = 
(lb — hb)/Ith- The SETOs regime exists only in the limit 
p > 1. 



One can readily verify that in the limit of an ideal cur- 
rent source the charge time dependence around each t n 
is Q(t) = Ib(t — ne/ It). The non-linear corrections in- 
stead add a stochastic dependence on the time evolution 
of the charge. The period, for instance, does not depend 
on (St) anymore for k — > oo, as is clear from Eq. (|18|) . 
This extreme limit is not realistic and thus the second 
term given in (|18[) is normally important. 

Reducing the current bias, the saturation charge Q s 
becomes of the order of c/2 (k = 2Q s /c — 1 < 1) and the 
non-linear behavior of Q(t) begins to correlate different 
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hopping events. The charge time evolution in this regime 
is characteristic and resembles a shark fin, as shown in 
Fig. (j4}. It is also clear from (fT8|) . that in this regime the 
stochastic fluctuations have the greatest impact on the 
average SETOs period. 

These two regimes can be identified on the current 
plot Fig. (IS]) and in the analytical expression (f2"Tj) as the 
two branches joined by a minimum of the voltage. The 
large bias behavior (k ^> 1) of (|2Tj) gives e y/im/p /(2C), 
which is the Averin-Likharev expression for the cur- 
rent (Ql ((/,/) ~ lb in this limit), while in the opposite 
limit the long exponential charging time is dominating: 
(^)«(e/C)^F/(y^ln(2/«)). _ 

The overall situation is summarized in Fig. (|7]l . Since it 
is very difficult in practice to experimentally reach large 
values of p. the plot in Fig. ([7]) suggests that a good 
experimental choice can be k = I, for which the SETOs 
appear for the lowest values of p. We will discuss in the 
following the correlation function of the charge in order 
to analyze quantitatively the evolution of the accuracy of 
the SETOs. 



IV. CHARGE-FLUCTUATION SPECTRUM 

A quantitative measure of the accuracy of the SETOs 
is given by the time correlation of the charge. This can 
be defined as 

S(t) = (Q(t + r)Q(f)) - (Q(t + t)) (Q(t)) , (22) 

where the average is performed over a statistical ensemble 
and the result does not depend on t, since the stochastic 
process is stationary. Note also that this quantity is pro- 
portional to the spectrum of current fluctuations through 
the shunt resistance R s : Sr{t) = S(t)/t^. In the case of 
voltage- biased junction (see Fig. ([T])), Sr gives the cur- 
rent fluctuations that can be directly measured through 
the load resistance R. 

For perfectly periodic charge oscillations the Fourier 
transform of ([22]) is given by a sum of Dirac delta func- 
tions at us = 2xn/T, with n integer. The non-periodic 
fluctuations introduce a finite width of these peaks. The 
form of S(ui) measures thus directly and with a simple 
procedure the accuracy of the periodic charge transfer. In 
this section we derive an analytical expression for S(oj) 
that allows us to better understand the origin of the fluc- 
tuations. In the next section we will compare these re- 
sults to those obtained numerically by Monte Carlo sim- 
ulations. 

In order to calculate the Fourier transform of S(t) 
for this stationary process it is convenient to define the 
charge Q(t) over a time < t < A, with A > T so that 
many SETOs are present in a single sample of Q(t). One 
can then calculate the Fourier series 

Q P = T Q(t)e v ^ , Q(*) = 5>-^*q p . 

Jo A p 

(23) 



Substituting Eq. (|23[) into Eq. (f22j) and averaging over t 
one obtains: 

S{r) = Y,W)e- ip2itT -Ql, (24) 
v 

which can be used numerically to compute the correlation 
function from the Monte Carlo data, or analytically, by 
performing the limit A — > oo. In particular the Fourier 
transform can be defined as 

,+A/2 

S{w)= / dTe WT ~ |r| S(r), (25) 

J-A/2 

which gives 

S(u) = Y,((\Qp\ 2 ) - Sp.oQD^Siu - Up ) , (26) 
p 

with ojp — 2irp/A. The presence of the Dirac delta func- 
tions is an artifact due to the periodic extension induced 
by the Fourier transform. In practice, since the frequency 
scale I /A is infinitesimal one can obtain the smooth func- 
tion S(u>) by averaging the expression (f2l))) for each value 
of us over a small interval 2n/ A. This simply gives that 

SH) = A(|Q P | 2 ) (27) 

for p^O (Wicncr-Khinchin theorem) . 

The problem is now reduced to the calculation of the 
Fourier series of the charge. Using the definitions of t n 
and T n given before in Eq. (jlip and assuming that the 
cxtrcma of the time evolution of Q(t) coincide with the 
two hopping events at times To and tjv we can write: 

^ / >T "+ i at 

Qp = Y, -Q{t)e^= (28) 

n=0 Jr " 

n=0 J o A 

In the limit of well-established SETOs the integral gives 
a contribution that fluctuates very little. On the con- 
trary the exponentials are much more sensitive to even 
small fluctuations of the tunneling times, since the phase 
results from the accumulation of many different hopping 
events. For this reason we expect that the upper inte- 
gration limit can be substituted with the period of the 
SETOs t„ + i — r n w T and we use Q/(e/2, t- (T- {St))) 
as the average charge dependence. The Fourier transform 
then takes the form: 

S{u) = NT{\F(u)\ 2 )A(uj), (29) 

where 

1 N ^ 

n=0 
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and 



A(tu) 



^Q f (e/2,t-T+(Sr))e iut 



(31) 



The quantity A can be readily evaluated: 



1 



(e %uT - 1) - 



-(Sr)/r B 



( e i"T _ e 7/r.j 



T/r s 



(32) 



In order to proceed we have to evaluate also the average 
of F(uj). It is convenient to express the time at which one 
event happens as a sum over the delays between previous 
events using Eq. (IT2l : 



F(u) 




(33) 



Now the average of |-F(u;)| 2 can be performed using the 
distribution function V(6t): 



(\F(u)\ 2 ) 



lfl + 2Re(^) 
N\ 1 1 - g(u>) 

where we introduce the quantities: 

g(u) = (e'^^A = 



SF_ 



d(ST)V(ST)e^ {Sl 



(34) 



(35) 



and 5F = 2Rc {g(g N — 1)/(1 — g) 2 }, whose contribution 
to S(u>) vanishes in the limit N — > oo. In conclusion we 
find: 



S(lo) 



l + 2Rc 



1 - g(u) 



TA(u) , 



(36) 



which constitutes the central result of this section. 

Wc are now in the position to study the spectrum of the 
charge fluctuations for the system at hand. As it can be 
seen from the form of (f34|) the function has a singularity 
for g — > 1. Since 

/•OO ft 

\g(uj)\ 2 = dt dt' V(t)P(t')2cos(ut) < 1 , (37) 



a singularity is present for uj — > or when the fluctuations 
are negligible so that g rj e l "(5( 5T )) . This picture predicts 
a series of peaks for the frequencies fl n = 2im/T with n 
integer. Small fluctuations introduce a finite width, regu- 
larizing the correlation function. The numerical integra- 
tion in the expression of g(u>) is straightforward, but it 
is also possible to obtain an analytical expression. Deep 
in the SETOs regime V(St) has a gaussian behavior and 
provides a short cutoff time, so that we can expand the 
exponential in (|35p to second order in St. This gives 



e fart. x 



(38) 



2iuj 



(St) 



n 2 iu + 4lo 2 t s (St 2 ) 
2(2 + k) 2 ~ 




FIG. 8. The relative width at half height of the first noise 
peak as a function of k, given by Eq. for p = 10. 



The maximum of Re{<7/(1 — g)} takes place for 
Arg[g(cj)] = that at lowest order in (St) gives for the 
position of the poles: 



U V (2 + «)i* 



(39) 



coinciding at linear order in (St) /t s with 2im/T. The 
phase of g (g = \g\e Z1 ^) thus vanishes at the minimum of 
Rc{g/(1 — g)}, so that near this point one can write at 
lowest order 4> ~ T(cj — The relevant part of (|3"4")l 

then reads: 



Re 



(l-|5l 



(l-\g\) 2 + 4> 21 



(40) 



and Ec|. (J29[) takes the simple Lorentzian form 



5(w) =i >i(n n ) 



r n /2 



Tl/A+(w-a n f 



(41) 



with the full width at half maximum r„ defined by: 



1-lfll 
T 



ml 



T{2 + K ) 2 



(St 2 ) (St)'- 



(42) 



The presence in this formula of the mean squared vari- 
ance of delay in the tunneling time, ( (St — (St)) 2 ), 
clearly indicates that the spread in the hopping events 
controls the width of the peak, as is physically expected. 
One also sees that the width of the poles increases with 
n. Performing the average with the gaussian distribution 
wc find ( (St - (St)) 2 ) /t 2 = 4(l-7r/4)/(/tp) and for the 
relative width T n /Q n we have the explicit expression: 



327rn(l -tt/4) 



In = . 

p K (2 + K ) 2 In 2 [(2 + k)J 



(43) 



Eq. (|4"3"|) allows to study the width of the peak in the 
charge-fluctuation spectrum, and thus, the accuracy of 
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the SETOs. It correctly gives that T 1 /(2tt/T) < 1 in 
the SETOs regime (1/p -C K -C p). It also shows (see 
Fig. ([8])) that the relative width T\/Vl\ is a monotonic 
decreasing function of the bias current (k) for a given 
value of the resistance (p) . 

From Eq. (|41|) it is clear that within this approxima- 
tion the weight of the Lorcntzian peak is controlled only 
by the form factor _4(f2„). Using its explicit expression 
one finds that for large p it becomes independent of p: 
A(n ± ) ~ e 2 / (4tt 2 + In 2 [2±*]). The full variation takes 
place in the K < 1 region, where the shape of Q(t) evolves 
from the shark-fin to the sawtooth form. From the form 
of Eq. ([43]) one can also see that for small k the peak 
broadens and the SETOs are washed out when Ti ~ f2i 
(for k ~ 1/p)- For large n the theory instead predicts 
that the relative width decreases monotonically [38| : 



T„ _ 87rn(l - tt/4) 
fl n pn 



(44) 



In this limit the SETOs disappear by a decrease of the 
weight of the peak, but within our approximation this 
is not seen. Actually for sufficiently large current bias 
(k ~ p) there is a finite probability that a single tunnel 
event is no more sufficient to bring the charge back in 
the Coulomb blockade region. This is quantified by the 
value of P n (t n ), which, contrary to our hypothesis, can 
become smaller than 1. One can expect the theory to 
roughly remain valid for the fraction of tunneling events 
that leaves Q(t n ) < e/2. This describes a peak that 
remains sharp, but that vanishes in weight as P n (t n ). 

Another interesting and relevant limit is the low- 
frequency behavior of S(u>). Expanding Eq. (|35p in ui 
one can show that 



N(\F(oj)\ 2 ) 



n 

2 AF1F2F3 



(45) 



F 2 F 4 



3F* 



12F- 4 



+ 



where F n = (J™) . If the fluctuations are negligible then 
F n = Fi and the noise at low frequency vanishes. In 
particular in our case this average takes a simple form if 
the explicit expression of $ is used: 



A(|^H| 2 ) 



(Sr 2 ) - (St)' 



1 



~L2~ 



(46) 

We thus find that in the SETOs regime the low-frequency 
noise is suppressed. Eq. (|4l)|) together with the ex- 
pansion of Eq. (|32[) for u> — > allows us to evaluate 
the zero- frequency Fano factor: F = Sj(0)/e{lj) = 
S , (0)/r 2 e (Ij), where Sj(uj) is the noise spectrum of the 
current through the junction and the relation S(uj) = 
5j(w)r 2 /(l +o; 2 r 2 ) holds exactly. The reduction of the 
current fluctuations in the large-p limit naturally leads 
to sub-Poissonian noise (F < 1) vanishing for p — > 00: 



(4-TT)(-2 + (K+l)U/T s f 

P k(2 + kY(U/t s Y 



(47) 



Nevertheless this expression is only qualitatively correct, 
since the analytical theory has been designed to describe 
accurately the noise for frequencies around the peak of 
the SETOs. This will be shown in the next section, where 
we present numerical simulations. 



V. NUMERICAL SIMULATIONS 

In this section we show numerical results obtained by 
Monte Carlo simulations for the charge-fluctuation spec- 
trum. The purpose of this section is to compare these 
results with the analytical calculations of the preceding 
sections, valid for p > 1 and 1/p -C k -C p, and to ex- 
plore the crossover region where the SETOs disappear. 

The Monte Carlo simulations are performed by gener- 
ating different realizations of the stochastic time evolu- 
tion of the charge Q(t) over a time much longer than the 
SETOs period. The time evolution of the charge is ob- 
tained by discretizing the time on a nonuniform grid, such 
that in the time interval At the charge varies by a small 
quantity and P = T(Q(t))At < Af < 1. The tunnel- 
ing event is accepted or refused by generating a random 
number between and 1 and by comparing it to P. The 
deterministic evolution of Q(t) between two events is sim- 
ply given by Eq. (j3)). The sequence of time intervals of 
deterministic evolution interspersed by tunneling times 
so constructed gives the full knowledge of Q(t) and con- 
stitutes the stochastic run. The square modulus of the 
Fourier transform of the charge is then easily calculated 
analytically piecewise, interval by interval. To obtain 
the noise as from (|27p . just a further average over several 
runs is needed. Typically J\f = .01, each run counts 10 3 
tunneling events, and an average over 10 4 realizations is 
performed. 

Let us now discuss the numerical results. We begin by 
comparing the form of the first peak in the noise spec- 
trum. It is shown scaled by the analytically calculated 
width Ti in Fig. ([9]) for different values of p and given 
k. The agreement is excellent. We then compare the full 
oj dependence of S(uS) obtained from Eq. (|36[) with the 
one calculated numerically. We show the comparison in 
Fig. dTDJ) for the relevant case k = 1 (which corresponds 
to the widest extension in p of the SETOs region) and for 
p ranging from 200 to 10. As expected, the agreement is 
very good deep in the SETOs regime, for p > 1, and at 
p = 10 all the essential features of the peak at ui = fli are 
still fairly well represented by the theory. The analytical 
approximation correctly finds the main contribution to 
the noise, but it tends to underestimate it far from the 
peak maximum. 

Let us now investigate how the SETOs disappear. 
From the experimental point of view a simple parame- 
ter that can be varied continuously is K ~ We thus 
plot in Fig. (TTT|) the evolution of S(lj) for given p = 10 
and 100 as a function of k. These two plots show several 
interesting features. The first striking one is the reduc- 
tion of the relative widths of the peaks by increasing k. 
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FIG. 9. (Color online) The noise peak at u — ttj obtained 
from Monte Carlo numerical simulations for values of p rang- 
ing from 80 to 200 in steps of 10. On the z-axis the fre- 
quency is shifted by Qi and scaled by the width Fi as given 
by Eq. (|42|) . on the y-axis the noise is scaled by T p. 




l 1.5 



FIG. 10. (Color online) Comparison between Monte Carlo 
data for the charge-fluctuation spectrum and the analytical 
results obtained from Eq. (|36p for different values of p, at 
fixed bias condition k = 1. 



This is predicted by the analytical expression (j43j) and 
the numerical calculations assess its validity even outside 
the region of applicability of the analytical theory. Note 
that in Fig. (|TTj) the frequency axis is rescaled with fii, 
thus the apparent weight of the peaks is reduced by the 
scaling, but it saturates in the region k <C p as predicted 
by the analytical theory, and then starts to decrease in 
the crossover region. The second visible feature is the ap- 
pearance of a wide Lorentzian zero frequency peak that 
remains the only structure for k 3> p. This structure 
is due to the charge noise induced at the capacitance 
by the Poissonian current fluctuations generated by the 
tunnel junction. By solving the electromagnetic prob- 
lem one finds: S(lu) = C 2 \Z(uj)\ 2 e (Ij), where e(I,j) is 
the standard tunnel junction Poissonian white noise and 
Z(uj) = R\\/{l + iLuR\\C), with R\\ = R s R t /(R s +R t ), the 




FIG. 11. (Color online) Monte Carlo spectra for different 
bias conditions k at fixed junction environment: p — 100 
(p = 10) in the upper (lower) graph. Increasing k means here 
moving along a horizontal line in Fig. toward the high- 
bias boundary of the SETOs region and it allows us to see 
how SETOs disappear. 



impedance between the current source and the voltage at 
the capacitance. This gives: 



S(u>) 



e(Ij)C 2 R 2 s 



l + 2p + p 2 (l + lo 2 R 2 C 2 ) 



(48) 



which fits our data very well (not shown) . 

Finally we show in Fig. (fTi?j) the evolution of S(u>) for 
k = 1 and for p evolving from 1 to 100. This figure gives 
an idea of the expected spectrum at the optimal value 
k = 1 as a function of p. It turns out that already at 
p = 3, S(uj) presents a very broad maximum and p = 5 
is probably sufficient to observe a clear structure in S(uj). 



VI. CONCLUSIONS 

In this paper we have studied theoretically electronic 
transport in a tunnel junction in the presence of a large- 
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FIG. 12. (Color online) Monte Carlo spectra for fixed bias 
conditions k = 1 and different values of p: The crossover 
from the SETOs regime to the p < 1 region is shown. 



resistive environment. The phenomenon of SETOs has 
been predicted to appear in this system for essentially 
infinite value of the external resistance, so that the junc- 
tion can be seen as current biased. We investigated under 
which conditions the SETOs appear for a realistic finite 
value of the environment resistance. We found analytical 
expressions for the current [Eq. (f20|) ] and for the charge- 
fluctuation spectrum [Eq. (J^U) with Eq. (|42[) ]. Our ana- 
lytical results describe very well the form of the peak in 
the charge noise, which can be regarded as the hallmark 
of the SETOs, since it quantifies the accuracy of the pe- 
riodicity in the charge time dependence. We find that a 
ratio of R s /Rt as low as 5 can be sufficient to observe a 
clear structure in S(lo) if the bias current is chosen such 
that h ~ 27 t h = e/R s C (this can be converted on a 
condition on the voltage bias V& ~ c/C). This ratio can 
be obtained experimentally and thus in principle SETOs 
can be observed through the measurement of the current 
or charge noise. 

The conclusion on the possibility of observing the ef- 
fect is thus optimistic. One should however keep in mind 
that, following the literature on this problem, the the- 
ory presented holds at low temperature (fc^T <C e 2 /2C) 
and neglects quantum fluctuations of the electromagnetic 
modes of the environment. Thermal and quantum fluc- 
tuations are expected to have a non-negligible influence 
on the transport mechanism at play in the system. Their 
evaluation requires a different technical approach and is 
beyond the scope of the present paper. 
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Appendix A: CALCULATION OF THE I-V 
CHARACTERISTICS IN THE SETOs REGIME 
WITH THE MASTER EQUATION APPROACH 

In this Appendix we find for the I-V characteristics 
expressions which take into account the spread in the 
distribution probability of the charge. A full statistical 
description of the behavior of the system can be given 
in terms of the probability a n {Q,t) that at time t the 
charge in the capacitance is Q and n charges have crossed 
the junction. Conservation of the probability and the 
master equation describing electron tunneling lead to the 
following set of coupled partial differential equations: 



d<r n (Q,t) = _3_ 
dt dQ 



(Q-Qs 



-a n (Q,t) 



(Al) 



-T(Q/C)a n {Q, t) + T((Q + e)/C)<r n _i(Q + e, t) . 

Eq. (1A1|) is a generalization of the equations given in 
Ref. [7| to include the information on the number of elec- 
trons which have tunneled. The general solution of this 
equation is difficult in the presence of a finite resistance. 
But in the SETOs regime we can find an explicit solu- 
tion by exploiting the fact that at every cycle the charge 
passes through the blocked range (— e/2 < Q < c/2). Let 
us assume that at t = the distribution function is: 



a n {Q,0) = 5(Q-e/2)5 n , no 



(A2) 



The differential equation (|A1|) for <r na (Q,t) can be easily 
solved, since it decouples from the others (cr no _i = 0): 



<r no (Q,t) = P(t)s(Q-Q f (^,t 



(A3) 



where P(t) = P n (i) as given by Eq. (|16p with t n = 
and Q / is defined in Eq. @ . Once we know the solution 
for <7„ we can substitute it into Eq. (|A1[) and find the 
solution for fj no+ \. This can be done by using the ansatz: 

°n 0+ i{Q, t) = J dQ'S (Q - Q f (Q', i)) f(Q', t) , (A4) 
that gives 

<rr»+i(Q,t) = f(z(Q),t)e* , (A5) 

with z{Q) = (Q - Q s )e^ + Q s , and f(z(Q),t) = for 
z{Q) < e/2 and t > t s In (l - . The resulting dif- 

ferential equation for / reads: 

= (Q, + I + (=«) - Q.)e-*) x (A6) 
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The equation can be integrated and gives for a no +i(Q, t) 

a no+1 (Q,t) =e^ P (r 3 In Q - *®.y) x (A7) 

r s (|-Q a )(z(Q) + |) 
ei? t C (z(Q)-§) 2 

that is non vanishing for Qi(t) < Q < Q2(t) where 
Qi(t) = Q s -e - (Q s - e/2)e"*/^ and Q 2 (t) = Q s - 
(Qs + c/2)e _ */ Ts . It is interesting to note that the dis- 
tribution has now a finite spread in Q induced by the 
combined action of the stochastic fluctuations and of the 
finite value of the resistance. This is in contrast with the 



simpler approximation used in the text to evaluate the 
average in Eq. (f2T)|). where we assumed that the spread 
was negligible, and that a delta function could be used 
to describe the distribution cr, i0 +i. 

We can now calculate the average charge on the junc- 
tion during a single oscillation: 

r T dt ( /c \ r Q2{t) \ 

(Q) = J j\Q f [-,t)P(t) + J dQQa no+1 (Q 7 t)\ . 

(A8) 

The numerical integration of this expression leads to the 
result shown in Fig. ([5]). 
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